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Abstract 



The ordinary Levy motion is a random process whose stationary independent increments are 
statistically self - affine and distributed with a stable probability law characterized by the Levy 
index a, 0< a <2. The divergence of statistical moments of the order q > a leads to an 
important role of the finite sample effects. The objective of this paper is to study the influence of 

these effects on the self - affine properties of the ordinary Levy motion, namely, on the "1/a 
laws" , that is, time dependence of the g-th order structure function and of the range. Analytical 

estimates and simulations of the finite sample effects clearly demonstrates three phenomena: 
spurious multi - affinity of the Levy motion, strong dependence of the structure function on the 
sample size at q > a, and pseudo - Gaussian behavior of the second - order structure function 
and of the normalized range. We discuss these phenomena in detail and propose the modified 
Hurst method for empirical rescaled range analysis. 

PACS number(s): 02.50.-r, 05.40. +j 
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I. INTRODUCTION. 



By Levy motions, one designates a class of random functions, which are a natural general- 
ization of the Brownian motion, and whose increments are stationary, statistically self-affine 
and stably distributed in the sense of P. Levy [1] . Two important subclasses are (i) a- stable 
processes, or the ordinary Levy motion (oLm), which generalizes the ordinary Brownian mo- 
tion, or the Wiener process, and whose increments are independent, and (ii) the fractional 
Levy motion, which generalizes the fractional Brownian motion and has an infinite span of 
interdependence. 

The theory of processes with independent increments was developed beginning from the 
Bachelier's paper [2] concerning Brownian motion. However, the rigorous construction of 
this process and studies of properties of its trajectories were undertaken by Wiener [3]. 
The modern presentation of the general theory of processes with independent increments 
is contained in [4]. The theory of the processes with independent increments possessing 
stable distributions has begun its history from the already cited work [1] and, later on, was 
developed by other prominent mathematicians. In particular, the properties of extremes 
of a- stable symmetric processes were studied in [5]. The geometric properties of their 
trajectories were considered in [6]. The monograph [7] contains a modern presentation of 
the theory of a- stable processes. 

The Levy random processes play an important role in different areas of application for 
at least two reasons. 

The first one is that the Levy motion can be considered as a generalization of the Brow- 
nian motion. Indeed, the mathematical foundation of the generalization are remarkable 
properties of stable probability laws. From the limit theorems point of view, the stable 
distributions are a generalization of widely used Gaussian distribution. Namely, stable 
distributions are the limit ones for the distributions of (properly normalized) sums of inde- 
pendent identically distributed (i.i.d.) random variables [8]. Therefore, these distributions 
(like the Gaussian one) occur, when the evolution of a physical system or the result of an 
experiment are determined by the sum of a large number of identical independent random 
factors. An important distinction of stable probability densities is the power law tails de- 
creasing as , a is the Levy index, < a <2. Hence, the distribution moments of the 
order q > a diverge. In particular, stably distributed variables possess a non-finite variance. 

The second reason for ubiquity of the Levy motions is their remarkable property of scale 

- invariance. From this point of view the Levy motions (like the Brownian ones) belong 
to the so - called fractal random processes. Indeed, the objects in nature rarely exhibit 
exact self - similarity ( like the Von Koch curve), or self - affinity On the contrary, these 
properties have to be understood in a probabilistic sense [9,10]. The random fractals are 
believed to be widely spread in nature. A coastline is a simple example of statistically self 

- similar object [9], as well as the spot of the Chernobyl contamination in the nearest zone 
[11]. On the contrary, the trace of the Brownian motion is statistically self - affine. Several 
numerical algorithms were developed in order to simulate fractional Brownian motion [10,12 

- 15]. They allow one to model many highly irregular natural objects, which can be viewed 
as random fractals [10]. The traces of the Levy motions are also statistically self - affine, 
therefore, one may expect that they are also suited for modeling and studies of natural 
random fractals. 
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The stable distributions and the Levy, or Levy - like, random processes are widely used 
in different areas, where the phenomena possessing scale invariance (in a probabilistic sense) 
are viewed or, at least, can be suspected, e.g., in economy [16 - 19], biology and physiology 
[20], turbulence [21] and chaotic dynamics [22], solid state physics [23], plasma physics [24], 
geophysics [25] etc. In this respect, the problems connected with experimental data process- 
ing are of great importance. It is also necessary to develop different numerical algorithms 
which allow one to simulate Levy motion with the given statistical properties. This, in turn, 
allows one to improve the methods aimed at analysis and interpretation of experimental 
data. Recently, three models of the oLm were proposed [26]. They can be interpreted as 
"difference schemes" to approximate the evolution equation for the distribution density of 
the ordinary Levy motion. In our paper we employ the different approximation, which is 
based on using Gnedenko limit theorem along with the inversion method for generating 
random variables. With the help of this approximation we study the consequences of self 

- affinity of the motion, namely, the time dependence of the structure functions and of the 
range. We show that the finiteness of the sample size plays an essential role when the conse- 
quences of self - affinity are considered. We demonstrate both analytically and numerically, 
that the finiteness of the sample size violates self - affinity, thus giving rise to spurious multi 

- affinity of the oLm. Further, the second order structure function and the normalized range 
(just these quantities are widely used in processing a huge number of various experimental 
data [12,27]) possess spurious, "pseudo - Gaussian" time - dependence. This circumstance 
allows one to suggest that at estimating the second order structure function and normalized 
range from experimental data the "Levy nature" of them can be easily masked. In order to 
avoid pseudo - Gaussianity when studying the range, we propose a modified Hurst method 
for rescaled range analysis. 

The paper is organized as follows. In Sec. 2 we discuss the property of self - affinity and 
its consequences. In Sec. 3 we propose a model for numerical simulation of the oLm. In 
Sec. 4 the effects of a finite sample size are studied. In Sec. 5 we present numerical results. 
Finally, the conclusions are exposed in Sec. 6. 



II. SELF - AFFINITY OF THE ORDINARY LEVY MOTION 

Let us proceed with the self - affine properties of the ordinary Levy motion denoted 
below as L a (t). In this paper we restrict ourselves by symmetric stable distributions. The 
characteristic function of the oLm increments is [4,7] 

p a , D (k, r) = (exp [ik {L a {t + r) - L a {t))\) = exp(- J D \k\ a r) (1) 

Here a is the Levy index, < a < 2, and D > is a scale parameter. The probability 
density 

/elk 
— exp(-ikx)p a)D (k,t) (2) 

— oo 

is expressed in terms of elementary functions in two cases: 

(i) ordinary Brownian motion, L a (t) = B(t) , which has a = 2 and the probability 
density 
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P2 - DM) = 75B! exp 



X 2 



ADt 



and 

(ii) ordinary Cauchy motion, which has a = 1 and the probability density 

/ x Dt/ir 

Pl,D {X,t) = 



D 2 t 2 + x 2 

At \x\ — > oo the probability densities of the oLm have power law tails, 

, , Dt 

p a ,D[X,t) OC 

\x\ 

The increments of the oLm are stationary in a narrow sense, 

L a (h + r) - L a (t 2 + r) = L a (ti) - L a (t 2 ), (3) 
and self-affme with the parameter H = 1/a, that is, for an arbitrary k > 

L a (i + r) - L Q (0 1 {k- 1 /" [L Q (i + rer) - L a (t)}} , (4) 

where = implies that the two random functions have the same distribution functions. 
The exponent 1/a, by analogy with the theory of fractional Brownian motion [28], is named 
the Hurst index of the ordinary Levy motion. 

We consider two corollaries of Eqs.(l) - (4), which, again by analogy with the definitions 
of Ref. [28], may be called "1/a laws" for the structure function and for the range. 

We first consider the structure function of the oLm. A "1/a law" for the structure 
function can be stated as follows: 

for all < a < 2 

S,(r, a) = (\LJt + r) - MOD = { q °* 9< ", (5) 

where 

V(/i; a) — J dx 2 \x 2 f J exp(—ixix 2 - \xi\ a ), (6) 

— oo — oo 

whereas for a = 2 (ordinary Brownian motion) 

S q (r; 2) = r^V(q; 2) (7) 

for an arbitrary q. 

Equations (5) - (7) have a direct physical consequence for description of an anomalous 
diffusion. Indeed, for the ordinary Brownian motion the characteristic displacement Ax(t) 
of a particle may be written in terms of the second order structure function as 
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Ax(r)=5 2 1/2 (r;2) = v / 2r 1 / 2 , (8) 

where the prefactor \/2 is simply V 1 ^ 2 (2;2) . One may note from Eqs.(6), (7) that, for 
the normal diffusion 

5V^;2)ocr 1 /2 (9) 

at any q and, thus any order of the structure function may serve as a measure of a normal 
diffusion rate: 

Ax(r) ^ S l q /q (q;2) oct 1 / 2 , (10) 

if one is interested in time - dependence of the characteristic displacement, but not in 
the value of the prefactor. We remind that usually just the time - dependence, but not the 
prefactor, serves as an indicator of normal or anomalous diffusion [23]. In analogy with Eqs. 
(9), (10) it follows from Eqs. (5), (6) that the quantity S^ q (r;a) at < a < 2 and any q 
< a can serve as a measure of an anomalous diffusion rate: 

Ax(r) « Sl /q (q; a) oc r 1/a , 0<q<a<2. (11) 

Here we have the case of a fast anomalous diffusion, or hyperdiffusion. 
The second corollary of Eqs. (1) - (4) is for the range of the oLm. A "1/ct law" for the 
range can be stated as follows: 

R(t) = sup L a (t) - inf L a (t) = r 1/a R(l). (12) 



0<t<T 



0<t<T 



For the ordinary Brownian motion r x l 2 R (r) has a distribution independent of r [29]. 

For the oLm, < a < 2, the statistical mean of the range, 

(R(t)) , behaves as r x l a at 1 < a < 2, and turn to infinity at < a < 1. 

Both "1/a law" are studied in numerical simulation in Sec. 5. However, at the end of this 
Section we point once more to an important property of the moments of stable distributions 
with the Levy index a < 2: the moments of the order q > a diverge. This property, in turn, 
manifests itself in divergence of the q - th order structure function at q > a , see Eq.(5), and 
of the mean of the range at a < 1. However, at experimental data processing both quantities 
are finite due to the finiteness of a sample size. In fact, one may expect that in the case of 
the Levy motion the finiteness of a sample size has a stronger influence on the results than 
in the case of the Brownian motion. Therefore, the estimates of such influence are needed. 
We study this problem in Sec. 4 before discussing results of numerical simulation. 



III. A SIMPLE WAY TO APPROXIMATE ORDINARY LEVY MOTION. 

The process of constructing approximation to the oLm can be divided into two steps. 

Step 1. At the first step we generate random sequence of i.i.d. random variables possess- 
ing stable probability law. These variables play the role of increments of the oLm having 
the Levy index , < a < 2. The value a = 2 corresponds to the ordinary Brownian motion, 
hence in this case the sequence of independent increments is generated with the use of a 
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standard Gaussian generator. Since in Ref. [12] the sequence generated by the Gaussian 
generator is called " approximate discrete - time white Gaussian noise" , we call the sequence 
generated at < a < 2 " approximate discrete - time white Levy noise" . 

We generate approximate white Levy noise possessing characteristic function 

Pa,D(k) =< exp(ikx) >= exp(— D \k\ a ). (13) 

At < a < 2 they have power law asymptotic tails [8], 

/ \ „r(l + a) sin(7ro;/2) , 
p a , D (x) oc I 1 \ awioo. (14) 

7T \X\ 

In the literature there exist different algorithms for generating random variables dis- 
tributed with stable probability law. We only mention two recently proposed schemes, 
which use the combinations of random number generators [30] and the family of chaotic dy- 
namical systems with broad probability distributions [31], respectively. However, we believe 
that the ways of generating stably distributed variables and, then the Levy motion are not 
exhausted, and various simulation models are needed, each of them may appear to be useful 
when studying some particular problem. In this paper we propose a simple approximation 
based on the Gnedenko limit theorem along with the method of inversion. Indeed, among 
the methods of generating random sequence with the given probability law F(x) the method 
of inversion seems most simple and effective [32]. However, it is well-known fact that its 
effectiveness is limited by the laws possessing analytic expressions for hence the direct 
application of the method of inversion to the stable law is not expedient. In this connection, 
it is natural to exploit an important property of stable distributions. Namely, such distribu- 
tions are limiting for those of properly normalized sums of i.i.d. random variables [8]. To be 
more concrete, we generate the needed random sequence in two steps. At the first one we 
generate an "auxiliary" sequence of i.i.d. random variables {£.,•} , whose distribution density 
F'[x) possesses asymptotics having the same power law dependence as the stable density 
with the Levy index has, see Eq.(14). However, contrary to the stable law, the function 
F(x) is chosen as simple as possible in order to get analytic form of F _1 . For example, 

F{x) = 2idw-y x<0 ' (15) 



F(x) = 1 - — r, X > 0. 

v ; 2(1 +x a Y 

Then, the normalized sum is estimated, 



am}l c . , 
where 

/ \ V" 

/ 71 \ 



i m 

X = —n Z T t ^ (16) 



2r(a)sin(7ra/2) 
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According to the Gnedenko theorem on the normal attraction basin of the stable law [8] , 
the distribution of the sum (16) is then converges to the stable law with the characteristic 
function (13) and D — 1. It is reasonable to generate random variables having stable 
distribution with the unit D, with a consequent rescaling, if necessary. Repeating N times 
the above procedure, we get a sequence of i.i.d. random variables {X(t)} , t — 1, 2, N . 
This is an approximate discrete - time white Levy noise. 

In the top of Fig.l the probability densities p(x) for the members of the sequence {X(t)} 
(m = 30 in Eq.(16)) are depicted by black points for (a) a = 1.0, and (b) a = 1.5. The 
functions p a ,i(x) obtained with the inverse Fourier transform, see Eq.(13), are shown by solid 
lines. In the bottom of Fig.l the black points depict asymptotics of the same probability 
densities in log- log scale. The solid lines show the asymptotics given by Eq.(14). The 
examples presented demonstrate a good agreement between the probability densities for 
the sequences {X(t)} obtained with the use of the numerical algorithm proposed and the 
densities of the stable laws. 

We would like to stress that a certain merit of the proposed model is its simplicity. It 
is entirely based on classical formulation of one of the limit theorems and can be easily 
generalized for the case of asymmetric stable distributions. It is also allows one, after some 
modifications, to speed up the convergence to the stable law. These problems, as well as the 
comparison with the other algorithms seems to be the subject of a separate paper. 

Step 2. With using approximate discrete - time white Levy noise X(t) the approximation 
to the oLm is defined by 

L a (t) = £,X(T) (17) 

T = l 

In Fig. 2 the approximate white Levy noises obtained with the numerical algorithm pro- 
posed are depicted by thin lines at 4 different Levy indexes. The thick lines depict the 
sample paths, or the trajectories, of the approximation to the oLm. It is clearly seen that 
with the Levy index decreasing, the amplitude of the noise increases. The large outliers lead 
to large "jumps" (often named as "Levy flights") on the trajectory. 

IV. EFFECTS OF LIMITED SAMPLE SIZE. 

In this Section we discuss what is the time dependence of the q - th order structure 
function and the range, which are estimated from experimental data. 
We first give an estimate for the mode of the maximum value. 

Suppose we have a sequence {X(t)}, t = 1,2,. ..,N, of i.i.d. random variables possessing 

X 

stable probability density p a ^ (x) and cumulative probability P a ,i(X < x) — J du p a ,i(u) 

— oo 

. Then, P^, is the probability that all N terms of the sequence are less than x. This, in 
turn, implies that P^i is the probability that the maximum value of N terms is less than x. 
Therefore, 

Vn{x) = «!(*))' = NP^\x)p^{x) (18) 
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is the probability density of the maximum value in the sample consisting of N terms. 
Let X ma _ x (N) be the mode of maximum value, that is, the most probable maximum value. 
It obeys an equation 

Since p a ,i(x) oc x~ l ~ a for large x, we may give the following estimate for the mode: 

X max (iV) oc N 1 ^ , < a < 2. (19) 

With the help of Eq.(19) we are able to roughly estimate the diverging statistical moment 
of the sequence {X(t)}, t = 1, 2, N, as 

X max (AT) 

J dXX*p aiX {X)<xX££<xN*l°-\ q>a. (20) 
o 

Now we turn to the structure function. It can be written as 

oo 

S q (r) = (\AL a (r)\ q ) = J dAL a | AL a \ q p aA {AL a , r), (21) 

— oo 

where AL a (t) = L a (t + r) — L a (t) , and the probability density is given by Eqs. (1), (2). 
We may introduce a stochastic variable £, such that AL a {r) = ^r 1 /" , and rewrite S q , 

oo 

S 9 = T*'° J d^'PaAO, (22) 

— oo 

where the probability density for a new variable is given by Eqs. (13), (14). 

To estimate the integral in Eq.(22), we use Eq.(20) with an important note that here N 
is equal to T/r, T is the total length of the sample. Therefore, the integral can be estimated 
as (T/r) 9 / " 1 , and 

S q oc r T q/a ~\ q>a. (23) 

Thus, the effects of a limited sample size manifest itself in Eq.(23), which replaces "theo- 
retical infinity" in Eq. (5) for q > a . A particular case is the second order structure function, 
for which we have "pseudo - Gaussian" relation (see Eq.(8)), 

S^ 2 oc t 1 / 2 [33]. The linear r- dependence of the q - th order structure function was 
found in [34]. 

Equation (23) points to violation of self - affinity of the oLm. Indeed, self - affinity 
implies that the - exponent of the structure function depends linearly on q, see Eqs. (4), 
(5). On the contrary, the rough estimate of the finite sample effects demonstrates that the 
exponent does not depend on gat q > a . Thus, in experimental data processing one may 
expect that the exponent smoothly changes its slope, thus giving rise to a convex curve. 
Such convexity seems as an indication of multi - affinity, see Ref. [35]. However, in our 
case this behavior of the exponent is stipulated not by "intrinsic" reason, but instead by 
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the influence of a finite size of a sample of the self - affine process. That is why we call this 
effect " spurious multi - affinity" . 

Now we consider the finite size effects for the range of the oLm. In the empirical rescaled 
range analysis, that is, at experimental data processing or in numerical simulation the range 
of the random process is divided by the standard deviation (that is, the square root of the 
second moment) for the sequence of increments, 



1/2 

^2 = i^wf 



IT N 

1 ,„,. n 2 



J t=l 



after subtraction of a linear trend. This procedure, called the Hurst method, or the 
method of normalized range, "smoothes" the fluctuations of the range on different segments 
of time series, and is used in a great variety of applications [27]. H. E. Hurst was the first [36] 
who has collected large statistical material relating to water levels and other phenomena, 
which indicates that the observed normalized range do not increase (as it is expected for the 
ordinary Brownian motion) like the square root of the observational period r, but instead 
like a higher power. However, the Hurst method is not satisfactory for the oLm because of 
the infinity of the theoretical value of standard deviation. What one should expect when 
estimating the normalized range of the oLm from the finite sample size ? Since (X 2 ) oc r 2 /" -1 
, see Eq. (20), then 

a 2 oc t 1 / 0-1 / 2 , whereas R(t) oc t 1 /" , and, thus, 

^ oc r 1 / 2 . (24) 

0"2 

Therefore, we conclude, that at the empirical rescaled range analysis the Hurst method 
gives "spurious", "pseudo - Gaussian" time dependence. 

In order to get the correct exponent 1/a, and, at the same time, to smooth fluctuations 
of the range, we propose to modify the Hurst method by exploiting the a - th root of the 
a- th moment instead of standard deviation, that is, 



I T \ 1 / a 



EW)! Q • (25) 



J t=i 




Since it has only weak logarithmic divergence with the number of terms in the sum 
increasing, then one has 

(26) 

where H pa 1/a is the Hurst index for the oLm with the Levy index , and the bar denotes 
averaging over the number of segments (having the length ) of the sample path. 

The expediency of using modified Hurst method when studying the Levy motion can be 
explained as follows. In general case, the Hurst exponent in Eq.(4) contains information not 
only on the Levy index of the increment distribution, but also on long - time correlations 
between the increments. In case of the ordinary Levy motion the correlations between non 
- overlapping increments are absent, and the Hurst index is equal 1/a. If correlations exist, 
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the Hurst index for the Levy motion differ from 1/a, and this circumstance leads to violation 
of the "1/a laws" for the structure function and for the range. Therefore, when treating 
experimental data, it seems expedient to estimate with the use of increment distribution 
(there are different methods for estimating parameters of stable distributions, see, e.g., 
[37]) and, then by testing laws" for the structure function and for the range, to get 

information about the presence (or absence) of long - time correlations. 

In the next Section we verify numerically the self - affine properties and finite sample 
effects with the help of approximation described in Sec. 3. 

V. SELF - AFFINE PROPERTIES AND FINITE SAMPLE EFFECTS IN 

NUMERICAL SIMULATION 

We study numerically r- and T - dependence of the q - th order structure function, 

S q oc T ^ q) T a{q \ (27) 

where, as before, is the time argument of the structure function, and T is the sample 
size. According to the analytical estimates 

m - { f • :>: , • < 28 > 

whereas 

a(q) = S- - 1, (29) 
a 

see Eqs. (5) and (23). 

We first study //vs in the relation at q less than a. In Fig. 3 a typical example is depicted 
by crosses at fixed q = 1/2. The q/a curve is shown by primes. One can be convinced 
himself that the 1/a law is well confirmed at q smaller than the smallest Levy index in 
numerical simulation. Also in the figure /ivs a is depicted by black points at q — 2. It is 
shown that the second order structure function, being estimated from a finite sample, leads 
to the spurious "pseudo - Gaussian" behavior. At the clarifying inset we show S q vs rin a 
log - log scale with q and being fixed. The exponent // in the main figure is obtained for the 
fixed a as a slope of a straight line at the inset. 

In Fig. 4a, b we plot //vs q for the Levy index 1.2 and 1.7, respectively. The analytical 
estimate, see Eq. (28), is shown by solid line for q < a and by dotted line for q > a . The 
results of simulation are shown by black points. The arrows indicate the value q — a , at 
which the bend of theoretical curves occur. It is shown that the results of simulation are 
well fitted by the analytical curves. 

Figure 5 demonstrates S q versus sample size T in a log - log scale. The Levy index a is 
equal 1.2. The dotted lines indicate S q vs T at q — a (horizontal line), q = 2.5a and 4a 
(sloping lines). The black points indicate simulation results (4 points for each q). As in 
Fig. 4, we see a good agreement between experimental results and analytical estimates. It 
implies, that Eq. (23), seeming very rough estimate, nevertheless, accounts for the finite 
sample size effects for the structure function of the oLm. 
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Figure 6 demonstrates the application of the modified Hurst method to the sample path 
of the approximation to the oLm with the Levy index a — 1. In Fig. 6a fluctuations of the 
range (thin curve) and those of a a (thick curve) are shown for the case, when the total length 
of the sample is divided into 64 segments, each of r = 16 lengthwise. Below the variations of 
the ratio R/a a are depicted. It is shown that fluctuations of the ratio is much smaller than 
those of the span. This circumstance justify the use of the ratio in the empirical analysis. In 
Fig. 6b the rescaled span, see Eq. (26) is depicted versus the time interval by black points 
in a log - log scale. The slope of the solid line is equal to H = 0.9. 

In Fig. 7 the Hurst index obtained with using Eq.(26) is depicted by crosses, whereas 
the curve 1/a is indicated by primes. By comparing Fig. 7 with Fig.3 one can see that the 
"1/a law" is better fulfilled for the structure function of the simulated process than for the 
range. We also note that the same conclusion follows for "t h laws" for fractional Brownian 
motion [15]. A more detailed discussion of this problem require the use of the theory of 
extremes for the processes with stationary increments. This is beyond the scope of our 
paper. However, we point readers attention on the Hurst index which is obtained with the 
use of "traditional" ratio Rjoi and is shown by black points in Fig. 7. In accordance with 
the discussion of Sec. 4, it can be clearly seen, that the standard deviation, being used in 
empirical analysis of the oLm, "suppresses" the variations of the range, thus giving rise to 
the spurious value H m 0.5 for all Levy indexes. On the contrary, only "smoothes" the 
variations of the range, thus leading to the (nearly) correct value H — 1/a. 

VI. CONCLUSIONS. 

The ubiquity of the Levy motion rises the problems related to the experimental data 
processing. In particular, the fat tails of the stable distributions allows one to suggest that 
the effects of a finite sample size play an important role. To study these effects, we propose 
an approximation which is based on using Gnedenko limit theorem along with the method 
of inversion. It allows us to simulate the ordinary Levy motion and study the finite sample 
effects when estimating structure functions and the range. The results of simulation are in 
a quantitative agreement with theoretical estimates. In particular, we find, that the second 
order structure function being estimated from the ordinary Levy motion sample paths as 
well as the Hurst method both lead to spurious "pseudo - Gaussian" time dependencies. 
We also propose to modify the Hurst method of rescaled range analysis in order to avoid 
"pseudo - Gaussian" relation. 

At the end we note that introduction of correlations into the sequence of i.i.d. stably 
distributed random variables (with the help of the method used in Ref. [15] for the sequence 
of i.i.d. Gaussian variables) allows us to study the finite sample effects for fractional Levy 
motion. 
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FIGURE CAPTIONS 

Fig. 1. Probability densities (above) and their asymptotics (below) are indicated for the 
sequences of random variables generated with the use of the proposed numerical algorithm at 
the Levy indexes (a) a — 1, and (b) a = 1.5. The probability densities and the asymptotics 
of the stable laws are indicated by solid lines. 

Fig. 2. Stationary sequences (thin lines) and ordinary Levy motion trajectories (thick 
lines) at the different Levy indexes. 

Fig. 3. Plots of the exponent in Eq.(12) versus the Levy index a at q = 0.5 (crosses) and 
q = 2 (black points). The q/a curve is depicted by dashed line. At the inset the structure 
function S versus is shown in a log - log scale at a — 1, q — 0.5. 

Fig. 4. The exponent \x (see Eqs. (27), (28)) versus the order q of the structure function 
for (a) a = 1.2 and (b) a = 1.7. The analytical estimate, see Eq. (28), is shown by solid 
line for q < a and by dotted line for q > a . The results of simulation are shown by black 
points. The arrows indicate the value q — a , at which bend of theoretical curve occurs. 

Fig. 5. Structure function S versus sample size Tin a log - log scale for the Levy index 
a — 1 and different values of q: q = a (horizontal line), q = 2.5a and q = 4a (sloping lines). 
The black points indicate simulation results. 

Fig. 6. (a) Variations of the range (thin curve), of a a (thick curve) and of their ratio 
(below) at the different time intervals for the oLm with a — 1. (b) Rescaled range, see Eq. 
(26), versus time interval in log-log scale (black points). Solid line has a slope H = 0.9. 

Fig. 7. Plots of the Hurst exponent H vs a estimated with the use of Eq.(16) (crosses) 
and with the use of the traditional Hurst method (black points). The 1/a curve is depicted 
by dashed line. 
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